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ABSTRACT 


In  this  expository  paper  the  concepts  of  competing  and 
complementary  risks  are  defined  and  a  survey  of  recent  results 
in  the  area  is  presented.  Identif iability  of  distributions  , 
both  univariate  and  multivariate,  useful  in  reliability  and 
survival  analysis  is  considered. 
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1.  INTRODUCTION 

The  problem  of  identifiability  (or  of  nonidentif lability) 
arises  naturally  in  a  number  of  physical  situations.  The  problem, 
in  general  terms,  can  be  defined  as  follows.  Let  U  be  an  observ¬ 
able  random  variable  with  distribution  function  F  ,  and  let  FA 
belong  to  a  family  F  =  {FQ:  0eH}  of  distribution  functions  indexed 
by  a  parameter  6.  Here  0  could  be  scalar  or  vector  valued. 

We  shall  say  0  is  nonidentif iable  by  l)  if  there  is  at  least 

one  pair  (0,  0")  0  ^  Q' ,  0,  0'efi  such  that  F0(u)  =  FQ^(u)  for  all 

u.  In  the  contrary  case  we  shall  say  0  is  identifiable.  It  may 

happen  that  0  itself  is  nonidentif iable  but  a  function  y(0)  is 

identifiable  in  the  following  sense:  For  any  0,  neft,  F  (u)  =  F  (u) 

o  n 

for  all  u  implies  y ( 0 )  =  y(n)>  in  this  case  we  may  say  0  is  partially 
identifiable. 

Puri  (1979)  has  surveyed  some  examples  arising  in  the  general 
literature.  The  purpose  of  the  present  paper  is  to  present  a  survey 
of  results  available  in  the  area  of  competing  risks  and  complimen¬ 
tary  risks.  We  would  primarily  consider  the  cases  where  we  have  an 
underlying  parametric  model.  In  Sections  2  and  3  we  define  the  pro¬ 
blems  of  competing  and  complementary  risks  and  discuss  the  case  when 


2 


the  underlying  random  variables  are  independently  distributed.  The 
case  of  dependent  random  variables  is  considered  in  Section  4.  The 
problem  of  estimation  is  briefly  discussed  in  Section  5.  Finally, 
in  Section  6,  some  open  problems  and  other  related  areas  of  research 
are  pointed  out. 
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2.  INDEPENDENT  RANDOM  VARIABLES  WITH  IDENTIFIED  EXTREMUM 

The  problem  of  competing  risks,  in  its  simplest  form,  may 
be  described  as  follows:  Let  be  a  random  variable  with  distri¬ 
bution  function  F^x),  (i  =  l,2,...,p).  We  assume  that  X^'s 
are  not  observable,  but  U  =  min (X^ , . . . ,X  )  is.  We  would  like  to 
estimate  F^'s  given  the  observations  on  U. 

This  model  finds  interesting  applications  in  a  number  of 
fields  but  particularly  in  problems  of  survival  analysis  and 
reliability  theory.  Thus  in  problems  of  competing 

risks,  or  in  studying  reliability  of  complex  systems,  an  individual 
or  a  system  of  components  may  be  exposed  to  p  different  causes  of 
death  (failure)  where  X^  is  the  time  to  death  from  cause 

(i  =  1,2,... ,p) .  Although  one  would  like  to  know  about  the 
distribution  of  the  x^s  only  observations  on  U’s  will  be  avail¬ 
able.  Basu  and  Ghosh  (1980)  give  some  of  these  examples.  For 
other  examples  and  a  survey  of  the  area  see  Birnbaum  (1979)  and 
David  and  Moeschberger  (1978) .  An  interesting  ecological  applica¬ 
tion  is  given  by  Anderson  and  Burnham  (1976)  who  study  the  popula¬ 
tion  ecology  of  the  mallard  where  the  two  causes  of  death  are 
"hunting"  and  "natural  mortality." 

In  case  the  X^ 1 s  are  independent  and  identically  distributed, 
the  problem  is  a  trivial  one  since 

P  n 

Fytx)  =  n  Fi(x)  =  (Fx(x))p  (1) 

i=l 

where  F  =  1  -  F  is  the  survival  function  (reliability  function) , 
Fu(x)  is  the  cdf  of  U,  and  F^fx)  is  the  common  distribution  function 

of  the  x. ' s. 

i 
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In  the  general  case,  when  the  X^'s  are  not  identically  dis¬ 
tributed,  the  problem  of  identif iability  can  be  illustrated  using 
the  following  example. 

Example  1. 

Let  (i  =  1,2, 3,4).  be  independent  random  variables  and  let 
X.  ~  e(A.).  That  is  let  the  cdf  of  X.  be 

11  l 

-A  .  x 

P  (x)  =  1  -  e  1  ,  x  >  0,  A^  >  0. 

If  A^  +  A2  =  A^  +  A^,  then 

min(X^,X2)  and  minjX^X^) 
are  identically  distributed. 

It  is  possible  to  introduce  additional  random  variables  so 
that  the  enhanced  family  with  the  added  information  becomes  identi¬ 
fiable.  In  this  case  we  call  the  original  family  of  distributions 
"rectifiable . "  Let  I  be  an  integer  valued  random  variable 
(I  =  l,2,...,p).  (U ,  I)  is  called  an  identified  minimum  if  I  =  k 
when  U  =  min (X^ ,X2 , . . . ,X^)  =  .  In  the  absence  of  I ,  U  will  be 

called  a  non- identified  minimum. 

Basu  and  Ghosh  (1980)  introduce  a  dual  problem,  called  the 
problem  of  complementary  risks,  where  instead  of  observing  the 
minimum  (identified  or  non-identif ied)  one  observes  the  maximum 
V  =  max (X^ ,X2 , . . . ,X  ) .  Basu  and  Ghosh  (1980)  introduce  examples 
to  show  that  this  problem  also  occurs  naturally  in  survival  analysis 
and  reliability  theory.  If  the  X^ ' s  are  iid  then,  as  in  the  case 
of  minimum,  one  can  trivially  obtain  the  cdf  F  (x)  of  X  from  that 

A 


of  v. 
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We  now  consider  the  case  when  the  X^'s  are  independently  but 
not  identically  distributed.  The  joint  probability  distribution 
of  (U , I)  is  specified  by  monotonic  functions 


H,  <x)  =  P(U  £  x,  I  =  k )  ,  (k  =  1 , 2 , . . .  ,p)  . 


(2) 


Then  Berman  (1963)  has  obtained  the  following  theorem. 


Theorem  1.  The  set  of  functions  (H  (x) }  is  related  to  the  set 

—  Jv 

{Fj (x) }  by  the  functional  equations 


x 

r 


Hk(x)  = 


n  (1  -  F. (t) )dF.  (t) ,  k  =  1,2, ...,p. 
jA  3  * 


(3) 


The  solution  of  this  set  of  equations  is 


x 

f 


F,  (x)  =  1  -  exp{ 


,-1 


[1  -  £  H.  ( t)  ]  dH,  (t)  }, 
j  =  l  3  k 


(4) 


k  =  1,2,. . . ,p. 


x 

! 


Let  pk  =  P(I  =  k)  and  G^(x)  = 


gk(t)dt  be  the  distribution 


function  of  the  conditional  random  variable  =  U|I  =  k.  Hence 


Hk (x )  =  PkGk (x) • 

Nddas  (1970)  has  given  an  alternate  expression  for  Theorem  1  as 


and 


rkU) 


Pkgk(x)/[  ?  P. (l  -  G  (x))J 
K  K  j=l  3  3 


Fk (x)  =  1  -  exp(- 


r^ (t)dt) 


(5) 


(6) 
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That  is  F.  (x)  are  uniquely  determined  by  the  distribution  of  the 

fk(x) 

observable  pair  (U,I).  Where  r,  (x)  =  -  is  variously  called 

*  1  -  Fk(x) 

as  force  of  decrement  (or  of  mortality) ,  age  specific  death  rate, 
failure  rate  function,  intensity  function,  and  hazard  function. 

Since  the  minimum  U  and  the  maximum  V  are  related  through 
the  relation 


Max  (X1,X2, . . . ,Xp)  =  -  Min (-  -  X2,...,-  Xp) 

it  follows  that  similar  results  also  hold  for  the  maximum. 

The  next  natural  questions  are:  (a)  What  if  the  minimum 
(maximum)  is  not  identified?  Could  we  obtain  the  distribution  of 
Xk  from  that  of  U?  And  (b)  What  can  we  say  about  the  identifi- 
ability  of  the  Fk's  if  the  Xk 1 s  are  not  independently  distributed? 
We  will  discuss  these  next. 
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3.  INDEPENDENT  RANDOM  VARIABLES  WITH  NONIDENTIFIED  EXTREMUM 
In  case  the  extremum  U  is  not  identified,  one  can  still 
uniquely  determine  F^(x)  under  certain  conditions.  To  this  end 
Basu  and  Ghosh  (1980)  obtain  the  following  theorem. 

Theorem  2 .  Let  F  be  a  family  of  pdf  on  with  support  (a,b) 
which  are  continuous  and  are  positive  to  the  left  of  some  point  A 
and  such  that  if  f  and  g  are  any  two  distinct  members  of  F  then 

lim{ f (x) /g  (x) } 
x+a 


exists  and  equals  either  0  or  Let  X.,...,X  be  independent 

L  P 

random  variables  with  respective  pdf’s  f.,f 2,..., f  in  F  and 

j.  ^  p 

Y  , Y_,...,Y  be  independent  random  variables  with  respective  pdf's 

X.  £. 

belonging  to  F.  if  min{X.,...,X  }  and  min{ Y..  ,  .  .  .  ,Y  }  have  identi- 

I  p  i  q 

cal  distributions,  then  p  =  q  and  there  exists  a  permutation 

(k^ ,k2 , • • • ,^p}  (l,2,...,p}  such  that  the  pdf  of  Y^  if 

f,  (i  =  l,2,...,p).  Anderson  and  Ghurye  (1977)  proved  a  similar 
Ki 

theorem  for  the  maximum.  As  an  application  of  the  above  theorem, 
consider  the  following  examples. 

Example  2 .  F  is  the  family  of  normal  distribution 


4>  (x|  u,o) 


1 

✓iro 


a 


Now , 


4>  (x|  v2'a2> 

lim  - 

X-*~«  <P  (x|  M  ,0 ^ 


1, 

0, 


if 

U1 

=  U2,  ox  =  o2 

if 

°2 

<  or  a1  =  o2 

and 

y2  " 

M1 

if 

°2 

>  Oj_  or  o1  =  o 2 

and 

U2  " 

W1 

(7) 


00 
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Conditions  of  the  above  theorem  are  met.  Hence  the  distributions 
are  identifiable. 

Example  3 .  F  is  the  family  of  negative  exponential  distributions 

f^ (x)  =  Ae  Ax,  x  ?  0,  A  >  0 

Here 


f ,  (x) 

,  .  X2 
lim - 

x->-0  f.  (x) 

A1 


(8) 


Here  conditions  of  the  above  theorem  are  not  met  and  hence  the 
distributions  are  not  identifiable. 

Note,  however,  if  the  maximum  is  observed  both  normal  and 
exponential  distributions  are  identifiable. 

There  may,  however,  be  situations  when  the  conditions  of 
the  above  theorem  will  not  be  met  yet  the  underlying 
family  of  distributions  may  be  identifiable.  For  example,  Basu 
and  Ghosh  (1980)  have  shown  that  the  gamma  and  the  Weibull  distri¬ 
bution,  which  contain  the  exponential  distribution  as  special 
cases,  are  not  identifiable.  In  fact  they  prove  the  following 
theorems . 

Theorem  3.  Let  the  pdf  of  be  given  by 


f  ^  (x)  -  f  (x  ;oe  ,  ih  )  = 


-x/3 .  a . -1 
e _ 1  x  1 

pi'li  r(V 


a  . 

1 


0,  P. 


(i  =  1,2, 3,4; 


where  and  are  not  both  equal  to  one  and  and  are  not 
both  equal  to  one.  Let  X^  and  X^  be  independent  random  variables. 
Similarly,  let  X^  and  X^  be  independent.  If  the  distribution  of 
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min(X^,X2)  is  identical  with  that  of  min(X3,X4),  then  either 
(alfa2)  =  (a3#ct4)  and  =  ( 3 3 , 6 4 ) 

or, 

(al '  a2 )  =  ^a4,a3^  and  =  ^4'e3^’ 

Now  consider  the  case  of  the  Weibull  distributions. 

pi 

p.  p.  -x  /6 . 

Let  f . (x)  =  ^  x  1  e  1 ,  x  >  0,  (0. ,p.  >  0) . 

1  U  .  11 

1 

pi 

-x  /0.  _  p.  p.-l 

Here  F. (x)  =  1  -  F. (x)  =  e  ,f. (x)/F. (x)  =  g—  x  ,  and  as 

1  1  1  1  O  . 

x  ■+■  0 , 

fi(x) 

FTTx) 

3 

Theorem  4.  Let  Xi  ~  W(pi,0i),  (i  =  1,2 ,3, 4)  be  independent  Weibull 
random  variables.  If  the  distribution  of  min(X^,X2)  is  the  same 
as  that  of  min(X3»X4),  then  either 

(pl»8l)  =  '  ^2  '  ®  2  ^  =  ^p4 '  ®4^  ' 

or, 


<pi'V  - 


provided  p^  ?  P2 . 


(p4,04)  and  (p2'e2) 


(P3/63) 
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In  Theorem  4  we  exclude  the  case  p^  =  For  if 

=  p,  say,  the  problem,  after  using  the  transformation 
(i  =  1,2),  reduces  to  that  of  the  exponential  distri¬ 
bution  . 

Note  that  the  above  suggests  that  Theorem  2,  in  its  present 
form,  provides  a  sufficient  condition  for  identif iability .  How¬ 
ever,  it  is  not  a  necessary  one. 
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4.  DEPENDENT  RANDOM  VARIABLES 
We  next  consider  the  case  when  the  random  variables 
. ,.,X  )  are  dependent.  We  restrict  our  discussion  primarily  to 
the  case  when  p  =  2.  There  are  many  physical  situations  where  it 
is  desirable  to  test  the  assumption  so  frequently  made  that  the 
X^'s  are  independent.  It  is  therefore  natural  to  study  the  extent 
to  which  (U,I)  or  U  determines  the  joint  distribution  of  X’s  can 
be  identified.  To  this  end  Basu  and  Ghosh  (1978)  pointed  out  the 
difficulties  using  the  following  construction.  Let  F(x^,x2)  be 
the  joint  distribution  of  (X^X^  ,  Ffx-^x^  =  P(X^>x  ,  X2>x2)  ,  and 
FfU^x  )  =  3F  ( x^/ x2)  /3x^ ,  (i  =  1,2).  For  simplicity  assume  that 
the  density  f(x1#  x2)  >  0  for  all  (x^x^.  Let 


Gi(x) 


exp 


x  _  _ 

-/  -Fi (Z,Z) (F ( Z , Z ) )  aZ 


(10) 


and  assume  /  -F.  ( Z  ,  Z)  (F  (z  ,  Z )  )  "'"dZ  diverges  for  i  =  1,2.  Then 

—  OO 

G^(x)  =  l-G^(x)  is  a  distribution  finction  and  (U,I)  has  the  same 
distribution  whether  (X^X^  is  distributed  according  to  Ftx^x  ) 
or  G  (x^,x2)  =  G^(x^)  *G  (x2)  .  Thus  our  porblem  could  have  a  satis¬ 
factory  solution  only  if  F  is  known  to  be  a  well  specified  para¬ 
metric  family  of  distributions.  Similar  results  for  nonidentifi- 
ability,  in  the  absence  of  specific  parametric  models,  has  also 
been  considered  by  Miller  (1977),  Tsiatis  (1975),  and  Rose  (1973). 
Tsiatis  (1978)  further  illustrated  the  magnitude  of  this  problem 
with  some  actual  data. 
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No  general  result  on  identi fiability  for  dependent  random 
variables  is  currently  available.  However  Peterson,  Jr.  (1975) 
obtained  the  following  interesting  inequalities  when  p  =  2.  Let 


Q-L  (Xj) 

=  P 

<X1 

> 

X1 

and  X^  < 

X2> 

q2 (x2) 

=  P 

(x2 

> 

X2 

and  X2  < 

xi> 

Q^O) 

=  P 

(X1 

< 

x2) 

=  V 

Q2(0) 

=  P 

<X2 

< 

Xl) 

=  V 

Then 


(ID 


x2)  = 

+  Q2(X2) 

-  P (x1  < 

X1  <  x2  <  x2^ ’  if  <  x2 

= 

+  Q2(x2) 

-  P (x2  < 

X^  <  X1  <  x1)  ,  if  x1  >  x2  ' 

(12) 

Qx  [max(x1,x2)  ]  +  Q2  (max  (x1  ,x2)  ]  <  F(x1,x2)  <  Q^Xj)  +  Q2<x2), 
and 


Q1(X1)  +  °2(X1)  ~  Fl(xi)  '  Ql(xl5  +  p2 

Q 1 ( x2 )  +  Q2 ( x2 ^  s  F2^X2)  -  pi  +  Q2(X2) 


(13) 


3ivariate  normal  distribution 

In  case  of  specific  parametric  models  questions  on  identifi- 
ability  has  been  settled  on  a  number  of  distributions.  We  summarize 
some  of  these  resultrfv  Consider  first  the  case  of  bivariate  normal 
distribution.  Let  (X^,X2>  ~  BVN  ( iij  ,  ^ ,  o  2  ,  p  12)  .  Na'das  (1971) 

showed  that  if  0  -  p^2  <  1  the  distribution  of  the  identified 
minimum  of  a  normal  pair  determines  the  distribution  of  the  pair. 
Na'das*  proof  is  not  complete,  however.  Basu  and  Ghosh  (1978)  com¬ 
pleted  Na'das'  proof  and  extended  it  to  the  case  of  nonidentified 
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minimum.  Recently  Gilliland  and  Hannan  (1980)  gave  an  elegant 
solution  of  the  problem  for  the  general  case.  Their  results  are 
given  below. 


Let  (X^ •  X2)  ^BVNf  p  j,  ^2 » o ^  *0  ^ ,  p )  .  Assume,  without  any  loss  of 
generality,  that  0^  >  or  =  o ^  and  P^  -  P2*  Identification 
of  parameters  '  p2  '  °1 '  °2  '  anC*  p  w:*-th:*-n  this  restricted  family  is 
equivalent  to  identification  of  parameters  up  to  the  switch  of 


(p  ,0  )  and  (P2'°2)  i-n  unrestricted  family.  Consider  the 


following  subfamilies  of  bivariate  normal  distributions  defined 
by  the  following  additional  restrictions  on  the  parameters. 


V 


0  2  >  0  ,  |  p  |  <  1 


V  °1  >  °2  =  0 


N2:  °2  >  P  =  ~1 


(14) 


°1  >  °2  >  O'  p  =  1 


N, 


°1  =  °2  “  O'  P  =  1  or  undefined. 


Let  N  =  uN.  .  Let  F„  denote  the  distribution  of 

0  1  u 


U  =  mintX^.X^)  when  (X^jX^)  has  distribution  F. 


Theorem  5.  If  F,G  t  NQ  and  Fu  =  Gy,  then  F  =  G. 


Theorem  6.  Suppose  that  F , G  t  N  and  F„  =  G„.  Then  F,G  c  H.  for 

u  u  1 


some  i  =  0,...,4.  If  F, G  e  N  ,  i  =  0,...,3  then  F  =  G.  If 


F , G  e  N^,  then  F  and  G  have  the  same  values  p^,o^,  and  0 2  values 


but  arbitrary  p^  ?  p^  values. 
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Ho  result  for  the  general  case  is  available  as  yet. 

However  Basu  and  Ghosh  (1978)  showed  the  identif iability  of  the 
trivariate  normal  distribution  given  the  distribution  of  the 
identified  minimum  and, given  that  for  each  pair  of  random  variable 


X. 


1 


X.  with  correlation  cofficient  p. .  and  deviations 
D 


o  . 

i 


and  a  .  , 
3 


1  -  p. .  o./a.  >  0  (i,j 


1,2,3;  i  *  j) . 


Note  most  of  the  identification  problems  considered  so  far 
are  also  valid  when  the  maximum  V,  instead  of  the  minimum  U,  is 


observable.  For  if  X  = ^  I  ^ is  p  variate  normal,  so  is  -X  and 


maxfX^, . . . , X  )  =  -min (-X^, . . . ,-Xp) , 

and  thus  identi fi cation  problem  for  the  maximum  can  be  restated 
in  terms  of  the  identification  problem  for  the  minimum.  Also, 
any  bivariate  distribution  obtained  through  strict  monotone  trans¬ 
formation  of  normal  variables  will  be  identifiable.  The  bivariate 
lognormal  distribution  is  thus  identifiable. 

Next  we  consider  identi f iability  of  a  number  of 

bivariate  distributions  useful  in  reliability  theory  and  survival 
analysis.  We,  in  particular,  consider  several  bivariate  exponen¬ 
tial  distributions.  A  survey  of  some  of  these  distributions  is 
presented  in  Basu  and  Block  (1975) .  These  include  the  bivariate 
exponential  distributions  of  Marshall  and  Olkin  (1967),  Block  and 
Basu  (1974),  and  Gumbel  (1960).  Basu  and  Ghosh  (1978,  1980)  have 
considered  identif iability  of  these  distributions.  There  results 


are  summarized  below. 
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(a)  Marshall  and  Olkin  bivariate  exponential.  The 
tail  probability  of  the  Marshall-Olkin  bivariate  distribution  is 


given  by 

F(x1,x2)  =  exp[-A1x1  -  ^2x2  ~  X12  max(x1»x2^  * 

(15) 

xl'x2,Xl'X2  >  °'  X12  "  0  ’ 

=  0,  otherwise. 


I 


Here  all  parameters  are  identifiable  if  (U,I)  is  observed.  However, 
if  only  U  is  observed  the  parameters  are  not  identifiable. 

(b)  Block-3asu  model.  Here  the  joint  density  function  is 
given  by 


f(xlfx2)  = 


xxi<x2 


a12)/(ax 


x2> 


exp 


~X1X1 


-  (A. 


X12 ) 


X1  <  X2 


AA2(A1  +  A12)/(A1  +  a2) 


exp 


(X1  +  X12)xl  “  X2X2 


X2  "  X1 


where  A  =  \1  +  A2  +  A12>  Here  the  parameters  are  not  identifiable 
at  all.  Note  that  the  model  proposed  by  Freund  (1961)  is  also  not 
identifiable,  since  it  is  related  to  the  Block-Basu  model. 


(16) 


Because  of  the  underlying  physical  assumptions  neither 
Marshall-Olkin  nor  Block-Basu  is  considered  a  suitable  one  as  a 
physical  model  when  the  maximum  is  observed. 

(c)  Gumbel  model  I.  Gumbel  (1960)  proposed  two  bivariate 
exponential  distributions.  The  first  one  is  given  by 


F(x1,x2)  =  1  -  expf-Aj^)  -  exp(-A2x2)  +  exp(-A^x^  -  A2x2  -  A^2x^x2)  , 

X1'X2'A1'A2  >  0  '  A12  "  °* 

Here,  the  parameters  are  identifiable  if  (U,I)  is  observed.  However 
if  the  nonidenti f ied  minimum  U  is  observed  only  A^2  and  A^  +  A2 
are  identifiable.  If  the  nonidenti f ied  maximum  V  =  max(x^,x2>  is 
observable,  then  A^2  is  identifiable  and  is  identifiable 

up  to  a  permutation. 

(d)  Gumbel  model  II.  Here  the  distribution  function  is 
given  by 


F(x1fx2)  =  (1  -  exp  -A^Ml  -  exp  -  A2x2)(l  +  A-^exp  (-A-^x  -A  x2)  ) 


Here  if  U  is  observable,  A^2  is  identifiable  and  ( A ^ , A 2 )  is 
identifiable  up  to  a  permutation. 

(e)  3ivariate  Weibull  distribution.  Like  the  bivariate 
exponential  distribution,  the  bivariate  Weibull  distribution  can 
be  defined  in  a  number  of  ways.  Consider  the  following  bivariate 
Weibull  survival  function 


F(xlfx2)  =  P(X1  >  xx,  X2  >  x 2 ) 

P,  ?2  Pi  P2  (19) 

=  exp  -  A^  A  +  A2x2  +  A12  max(x1i,x2  ) 

Here  again  if  U  is  observable,  p^  and  p2  are  identifiable  up 
to  permutation.  Also  A^  and  A 2  +  A^2  or  A2  and  A^  +  A^2  are 
identifiable . 


r 
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5.  ESTIMATION  OF  PARAMETERS 

Estimation  of  parameters- based  on  (U,I),  the  identified  mini¬ 
mum  has  been  considered  extensively.  Basu  and  Ghosh  (1978)  have 
considered  the  problem  of  estimation  of  parameters  based  on  U  alone 
for  the  bivariate  normal  distribution.  Similar  results  can  be 
obtained  using  V. 

We  consider  some  other  models.  The  pdf  of  V,  assuming  inde¬ 
pendence,  is  given  by 

f  (t)  =  f.  (t)F.(t)  +  f_(t)F.(t).  (20) 

v  1  2  2  1 

The  paraim  t-ers  can  therefore  be  estimated  numerically  using  method 
of  mavxmum  xikelihood. 

As  an  illustration,  let  the  density  function  of  be  given 
by 


fi(xi)  =  ^  exp(-A^x^)  ,  x^A.^  >0,  (i  =  1,2). 

Then  the  density  function  of  V  is  given  by 

fv(t)  =  A^  exp(-A1t)  +  A2  exp  ( —  X  ^  t )  -  (A^  +  A2)  exp  -(A.^  +  A2)t. 


Hence  one  can  find  the  maximum  likelihood  estimates.  A  simpler 
method  of  estimation  is  by  the  method  of  moments.  We  can  readily 
show 


and 


(21) 


E(V2) 


2  +  2 _ 2 

Xj2  (Ax  +  X2)2 


(22) 


2 

Replacing  E(V)  and  E (V  )  by  their  simple  estimates  V^/n  and 
2 

/n  one  can  solfe  for  ^  and  2  by  iteration. 

Similarly,  suppose  X^,X2  are  dependent  and  follow  the  follow¬ 
ing  distribution  of  Gumbel. 


F (x^ ,x2 )  =  1  -  e 


-xixi 


-X2X2 


. -Xlxl-X2x2-X12xlx2 


Here  the  density  function  of  the  maximum  V  is  given  by 


-At  -A_t  -(A  +A  )t-A  t2 

fv(t)  =  Aje  1  +  A2e  2  -  (Aj.  +  A2  +  2A12t)e  1  '  (23) 

Using  this  one  can  again  estimate  the  parameters  by  the  method  of 
moments . 


< 
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6 .  CONCLUDING  REMARKS 


In  this  section  we  point  out  a  number  of  areas  in  which 
additional  work  is  being  carried  out  and  we  also  point  out  a  few 
open  problems. 

So  far  we  have  assumed,  in  the  competing  risks  problem,  that 

the  observed  random  variable  U  is  the  minimum  of  p  (unobserved) 

random  variables  X2,...,Xp.  In  analyzing  mortality  data  the 

above  is  interpreted  as  follows.  Consider  a  population  in  which 

K  causes  of  death,  C^ ,  C2 ,...,Cp,  are  operating.  Each  individual 

in  this  population  is  exposed  to  risk  of  dying  of  any  one  of  these 

causes.  One  can  recognize  two  kinds  of  distributions  associated 

with  death  due  to  cause  C. : 

1 

(a)  the  survival  distribution  F.  (t)  due  to  cause  C. , 

ia  i 

condi tionalxy  that  C^  is  the  cause  of  death,  in  the  presence  of 
other  causes; 

(b)  the  survival  distribution  F. (t)  due  to  C. ,  if  C.  is  act- 

i  l  i 

ing  alone. 


It  is  tacitly  assumed  that  corresponding  forces  of  mortality 
(failure  rates)  are  the  same.  Gail  (1975),  Elandt-Johnson  (1976), 
and  others  have  explored  the  implication  of  this  assumption. 

A  second  assumption  is  that  the  potential  survival  times  X^ ' s 
are  independently  distributed  with  continuous  distribution  function 
F^  (i  =  l,2,...,p).  Some  results  to  this  end  are  given  in  sections 
2,  3,  and  4.  For  another  interesting  direction  of  research  see 
Miller  (1977)  ,  Desu  and  Narula  (1977)  ,  Langberg,  Proschan,  and  Quinzi 
(1977)  (1978) ,  and  the  references  therein.  Desu  and  Narula  consider 

the  problem  of  estimating  the  distribution  function  F^(t)  =  P(X^  s  t) 


-r  ••  -  *>. 


i  ,  -*,?*.-**•-  ‘  •;  T" 
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and  provides  sufficient  condition  on  the  distribution  of 
(X^,...,X  )  under  which  such  an  estimation  is  possible. 

A  third  direction  of  research  has  been  the  interpretation 
of  competing  risks  problems  in  terms  of  some  stochastic  process. 

Chiang  (1968)  has  studied  the  problem  of  competing  risks  using 
time-nonhomogeneous  Markov  processes.  Clifford  (1977)  and  Berlin, 
Brodsky,  and  Clifford  (1977)  have  considered  the  problem  of  iden- 
tifiability  for  this  situation. 

The  problems  described  in  the  previous  chapters  can  be  ex¬ 
tended  in  several  directions.  The  author  is  currently  working  on 
some  of  the  problems  stated  below. 

(a)  As  mentioned  before,  Theorem  2  of  Section  3  does  not 
provide  a  necessary  and  sufficient  condition.  It  would  be  desir¬ 
able  to  improve  on  this  result. 

(b)  Most  of  the  identifications  results  obtained  so  fare  are 
for  the  case  of  two  competing  causes.  These  need  to  be  generalized 
for  the  case  of  any  number  of  variables.  Some  results  to  this  end  have 
been  obtained  by  Basu  and  Ghosh  (1980b) .  Necessary  algorithms  for  esti 
mating  the  parameters  should  also  be  obtained. 

(c)  The  concept  of  competing  risks  is  well  known.  Basu  and 
Ghosh  (1980)  first  coined  the  term  complementary  risks  for  the 
dual  problem.  In  reliability  theory  the  corresponding  problems  are 
for  series  and  parallel  systems  It  is  natural  to  pose  the  follow¬ 
ing  general  problem  corresponding  to  a  k-out-of  p  system  (k  <.  p)  . 


Recall  a  system  is  called  k-out-of-p  if  the  system  operates  so  long  as 
k  or  more  components  function.  Let  X(r)  be  rt^  ordered 
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Given  the  distribution  function  of  X^  could  we  uniquely  determine 
the  marginal  distribution  function  of  X.  (i  =  1,  2,...,p)?  Note, 
if  r  =  1,  we  obtain  the  case  of  competing  risks.  And  if  r  =  p,  we 
obtain  the  case  of  complimentary  risks.  Below  we  get  some  partial 
solution  for  the  case  of  three  independent  exponentially  distributed 
random  variables. 

Let  X^ ,  X2 ,  X^  be  three  independent  random  variables  and  let 
X^  ~  e(A^)  (i  =  1,  2,  3).  Suppose  distribution  of  X(2)  known 
and  it  is  known  which  of  the  three  variables  X^,  X^ ,  and  X^  is  X(2)’ 
That  is  we  assume  X,0.  to  be  identified.  If  the  X.'s  have  a  common 
distribution  F  with  density  function  f  we  can  again  solve  the 
problem  readily.  For 

F(2)  (t)  =  P(X(2)  >  t)  «  2(\  -  dpi  +  dpi)  (24) 

In  general,  if  X^  has  distribution  function  F^  and  density  function 

f .  we  have 
i 

f(2)(t)  =  (FirKff2(t)F3(t)  +  (F1(t)f2(t)F3(t)  (25) 

For  the  exponential  case  we  have  the  following  theorem. 

Theorem  7 .  Let  ,  X2 ,  X3  be  independent  random  variables  with 
X^  ~  e(A^)  (i  =  1,  2,  3).  Similarly  let  X^"*,  X2',  X^  be  independent 
random  variables  with  X/  ~  e(AT)  (i  =  1,  2,  3).  If  the  distribu¬ 
tion  of  the  identified  rfc^  ordered  statistics  (among  X^,  X2,  X3) 

and  (among  X',  X',  X')  are  identical  then  either 

(  A^  ,  A  2  ,  A  3 )  =  ( A  ^ ,  A  2 ,  A  3 )  or  ( A  ^ ,  A  2 ,  A3)  =  (A  3 ,  A  2 ,  A  ^ ) , 

Since  the  result  is  well  known  for  r  =  1  and  r  =  3  consider 
the  case  when  r  =  2 . 
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Proof.  Since  the  order  statistics  is  "identified"  let  us  assume 
=  X2  and  =  X' •  Now  the  densities  of  and  U"^  are 

identically  equal.  Hence 


-A_t  -At  -At  -(A-+A,)t 

,  2  r  1  ,  3  131 

\2e  Le  +  e  -  2e  J 


(26) 


-U0t  -M,t  -U,t  -  (U  +U  t  )  t_ 

=  u2e  2  [e  1  +  e  3  -  2e  1  3  J  for  all  t. 


or, 


A  2  {  (A  +\2)-(u1+u2)}t 

—  e 


.  ^  -  (A  -Xi't  -A  t 

,  ^  -rM3-ui>t  -u  t 

1  +  e  -  2e 


for  all  t.  (27) 

For  simplicity  assume  A  ^  <  A  ^  and  u  ^ .  Then  as  t  +  ®,  the 

i J. jht  hand  side  of  (27)  tends  to  1.  Hence  the  left  hand  side  limit 
of  (27)  is  also  1.  This  implies  u  =  X1  and  u2  =  A2>  Substituting 
these  values  back  in  (26)  we  note  u ^  =  A^. 

More  results  in  this  direction  have  been  obtained  by  Basu 


and  Ghosh  (1980b) . 
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